Dark matter, density perturbations and structure formation 



A. Del Popoloi'2.3 

^BogaziQi University, Physics Department, 80815 Bebek, Istanbul, Turkey 
^ Dipartimento di Matematica, Universita Statale di Bergamo, via dei Caniana, 2, 24-127, Bergamo, ITALY 
^Istanbul Technical University, Ayazaga Campus, Faculty of Science and Letters, 344^9 Maslak/ISTANBUL, Turkey 

Abstract — -This paper provides a review of the variants of dark matter which are thought to be 
fundamental components of the universe and their role in origin and evolution of structures and 
some new original results concerning improvements to the spherical collapse model. In particular, 
I show how the spherical collapse model is modified when we take into account dynamical friction 
and tidal torques. 



1. INTRODUCTION 

OOThe origin and evolution of large scale structure is today 
je outstanding problem in cosmology. This is the most 
idamental question we can ask about the universe whose 
^solution should help us to better understand problems as 
rae epoch of galaxy formation, the clustering in the galaxy 
distribution, the amplitude and form of anisotropics in the 
fjnicrowave background radiation. Several has been the ap- 
proaches and models trying to attack and solve this problem: 
r^toione has given a final answer. 

^e leading idea of all structure formation theories is that 
S^Hictures was born from small perturbations in the other- 
^se uniform distribution of matter in the early Universe, 
^rh ich is supposed to be, in great part, dark (matter not 
detectable through light emission). 

I ^W ith the term Dark Matter cosmologists indicate an 
hypothetic material component of the universe which does 
rrrrtt emit directly electromagnetic radiation (unless it decays 
particles having this property ([l|, but also see 0)). 
irk matter, cannot be revealed directly, but nevertheless 
sis necessary to postulate its existence in order to explain 
JlCi discrepancies between the observed dynamical proper- 
ties of galaxies and clusters of galaxies and the theoretical 
^^edictions based upon models of these objects assuming 
^gat the only matter present is the visible one. If in the 
f^ace were present a diffused material component having 
gravitational mass, but unable to emit electromagnetic 
. Eajdiation in significative quantity, this discrepancy could 
eliminated (0]). The study of Dark Matter has as its 
feadlity the explanation of formation of galaxies and in 
. funeral of cosmic structures. For this reason, in the last 
decades, the origin of cosmic structures has been "framed" 
in models in which Dark Matter constitutes the skeleton of 
cosmic structures and supply the most part of the mass of 
which the same is made. 

There are essentially two ways in which matter in the 
universe can be revealed: by means of radiation, by itself 
emitted, or by means of its gravitational interaction with 
baryonic matter which gives rise to cosmic structures. 
Electromagnetic radiation permits to reveal only baryonic 
matter. In the second case, we can only tell that we are in 
presence of matter that interacts by means of gravitation 
with the luminous mass in the universe. The original 
hypotheses on Dark Matter go back to measures performed 
by Oort ([3|) of the surface density of matter in the galactic 
disk, which was obtained through the study of the stars 
motion in direction orthogonal to the galactic plane. The 
result obtained by Oort, which was after him named 
"Oort Limit", gave a value of p = 0A5Mopc^^ for the 
mass density, and a mass, in the region studied, superior 
to that present in stars. Nowadays, we know that the 
quoted discrepancy is due to the presence of HI in the solar 



clusters (e.g. Coma Cluster) and the total mass contained 
in galaxies of the same clusters. These and other researches 
from the thirties to now, have confirmed that a great part 
of the mass in the universe does not emit radiation that 
can be directly observed. 



1.1 Determination of D, and Dark Matter 

The simplest cosmological model that describes, in a suf- 
ficient coherent manner, the evolution of the universe, from 
10~^s after the initial singularity to now, is the so called 
Standard Cosmological Model (or Hot Big Bang model). It is 
based upon the Friedmann-Robertson- Walker (FRW) met- 
ric, which is given by: 



ds^ = c^dt^ - a{tf 



dr^ 



1 — kr^ 



+ r'^{d0^ + sine'^dcj)'^) 



(1) 



where c is the light velocity, a(t) a function of time, or a 
scale factor called "expansion parameter" , t is the time co- 
ordinate, r, 9 and the comoving space coordinates. The 
evolution of the universe is described by the parameter a(t) 
and it is fundamentally connected to the value p of the av- 
erage density. 

The equations that describe the dynamics of the universe 
are the Friedmann's equations (Q) that we are going to in- 
troduce in a while. These equations can be obtained starting 
from the equations of the gravitational field of Einstein (Q): 



^ik ^ nQikTi — 7--'-ik 
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(2) 



where now, Rik is a symmetric tensor, also known as Ricci 
tensor, which describes the geometric properties of space- 
time, gik is the metric tensor, R is the scalar curvature, Tik 
is the energy-momentum tensor. 

These equations connect the properties of space-time to the 
mass-energy. In other terms they describe how space-time 
is modeled by mass. Combining Einstein equations to the 
FRW metric leads to the dynamic equations for the expan- 
sion parameter, a(t). These last are the Friedmann equa- 
tions: 



d{pa^) = -pd{a^) 



SttG 



2- + — + — = -SnGp 



(3) 



(4) 



(5) 



where p is the pressure of the fluid of which the universe 
is constituted, k is the curvature parameter and a(t) 
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the components of the today universe are galaxies. If we 
assume that galaxies motion satisfy Weyl ([^) postulate, 
the velocity vector of a galaxy is given by = (1,0,0,0), 
and then the system behaves as a system made of dust for 
which we have p = 0. Only two of the three Friedmann 
equations are independent, because the first connects 
density, p to the expansion parameter a(t). The character 
of the solutions of these equations depends on the value 
of the curvature parameter, fc, which is also determined 
by the initial conditions by means of Eq. 3. The solution 
to the equations now written shows that if p is larger 
than pc — 1^ — 1.88 * lQ~'^^g/c-m? (critical density, 
which can be obtained from Friedmann equations putting 
t = toi = 0, and H = 100km/ sMpc), space-time has 
a closed structure (fc = 1) and equations shows that the 
system go through a singularity in a finite time. This 
means that the universe has an expansion phase until it 
reaches a maximum expansion after which it recoUapse. 
If /9 < pc, the expansion never stops and the universe 
is open fc = — 1 (the universe has a structure similar to 
that of an hyperboloid, in the two-dimensional case). If 
finally, p — p^. the expansion is decelerated and has infinite 
duration in time, fc = 0, and the universe is flat (as a plane 
in the two-dimensional case). The concept discussed can 
be expressed using the parameter fl = In this case, the 
condition = 1 corresponds to k — 0, < 1 corresponds 
to fc = —1, and > 1 corresponds to fc — 1. ^ 
The value of O can be calculated in several ways. The most 
common methods are the dynamical methods, in which the 
effects of gravity are used, and kinematics methods sensible 
to the evolution of the scale factor and to the space-time 
geometry. The results obtained for fi with these different 
methods are summarized in the following. 

Dynamical methods: 

(a) Rotation curves: The contribution of spiral galaxies 
to the density in the universe is calculated by using their 
rotation curves and the third Kepler law. Using the last it 
is possible to obtain the mass of a spiral galaxy from the 
equation: 

M{r) = v^r/G (6) 

where v is the velocity of a test particle at a distance 
r from the center and M(r) is the mass internal to the 
circular orbit of the particle. In order to determine the 
mass M is necessary to have knowledge of the term in 
Eq. ^ and this can be done from the study of the rotation 
curves through the 21 cm line of HI. Rotation curves of 
galaxies are characterized by a peak reached at distances 
of some Kpcs and a behavior typically flat for the regions 
at distance larger than that of the peak. A peculiarity 
is that the expected Keplerian fall is not observed. This 
result is consistent with extended haloes containing masses 
till 10 times the galactic mass observed in the optical 
([To|V The previous result is obtained assuming that the 
halo mass obtained with this method is distributed in a 
spherical region so that we can use Eq. ([6]) and that we 
neglect the tidal interaction with the neighboring galaxies 
which tend to produce an expansion of the halo. After 
M and the luminosity of a series of elliptical galaxies is 
determined, the contribution to the density of the universe 



is given by p =< ^ > £ where £ is the luminosity per 
unit volume due to galaxies and can be obtained from 
the galactic luminosity function (j3{L)dL, which describes 
the number of galaxies per Mpc' and luminosity range 
L,L + dL. The value that is usually assumed for £ is 
£ = 2AhlO^LboMpc~^ . The arguments used lead to a value 
of rig for the luminous parts of spiral galaxies of Vlg < 0.01, 
while for haloes flh > 0.03 — 0.1. The result shows that 
the halo mass is noteworthy larger than the galactic mass 
observable in the optical f[ll|). 

(b) Virial theorem: 

In the case of non spiral galaxies and clusters, the mass 
can be obtained using the virial theorem 2T -1-^ = 0, with 

T - < vl > (7) 

where < > is the velocity dispersion along the line 
of sight. After getting the value of M of the cluster by 
means of the virial theorem one determines L by means of 
observations. Given M and L, the value of Q, for clusters 
is obtained similarly to the case of spiral galaxies. Usual 
values obtained for Vl are = 0.1 — 0.3 ([Hi). A problem of 
the quoted method is that in general the results obtained 
are the right one only for virialized, spherically symmetric 
clusters. In general, clusters are not virialized objects: even 
Coma clusters seems to have a central core constituted by 
more than one blob of mass ([l^). 

(c) Peculiar velocities: 

The velocity of a galaxy can be written as: 

Yg = Hr + (5v (8) 

where H is Hubble constant. The previous equation shows 
that the motion of a galaxy is constituted by two compo- 
nents: the velocity of the galaxy due to the Hubble flow 
and a peculiar velocity i5v, which describes the motion of 
the galaxy with respect to the background. In the linear 
regime, as we see in a while, we find that on average on a 
scale of length A, it is: 

(9) 

((l3|). Then given the overdensity ^ on scale A and 5v, 
it is possible to obtain il. The overdensity ^ can be ob- 
tained from the overdensity of galaxies using the rela- 
tion ^ = ^fe-i with 1 < 6 < 3. Using IRAS catalog in 
order to obtain the overdensity in galaxies one finds = 1. 
The values of obtained using the method of peculiar ve- 
locity assume that the peculiar velocity fields describe in 
an accurate way, the inhomogeneity in the distribution of 
underlying mass. We should note that the peculiar velocity 
method has some difficulties. In general, in order to obtain 
these last it is necessary to determine the redshift and the 
distance of galaxies and by using these data it is possible to 
obtain the peculiar velocity: 

Vpec = zc- Hod (10) 

It is evident that there are problems in measuring the 
distance d, problems connected to difhculties in finding 
trustable indicators of distance. Moreover the peculiar 
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(d) Kinematic methods: 

These methods are based upon the use of relations be- 
tween physical quantities dependent on cosmological param- 
eters. An example of those relations is the relation luminos- 
ity distance- redshift: 



(11) 



where is Hubble constant nowadays, z is the redshift, 
dh = the luminosity distance, L the absolute luminos- 
ity, and F the flux. By means of the relations luminosity- 
redshift, angle-redshift, number of objects-redshift, it 
is possible to determine the parameter of deceleration 
go = --fP^ (oo and Hq = = 100hkm./Mpcs are the 
scale factor and the Hubble constant, nowadays). At the 
same time go is connected to £7 by means of go = §j in a 
matter universe. One of the first test used, the luminosity 
distance-rcdshift has several problems due to effects of the 
evolution of sources. Uncertainties in the knowledge of the 
effects of galactic evolution on the intrinsic luminosity of the 
same has not, in the past, permitted to find deflnitive values 



of go. For this reason, introduced another kinematic 



test: number of galaxies- redshift. This test is based on the 
count of the number of galaxies in a comoving element of 
galaxies, defined by the surface dJl and the redshift dz. This 
number depends on go. Nevertheless the effects of evolution 
of sources influences on the results of the test, it is more 
sensible to the evolution of number of sources than to the 
evolution of luminosity, on which there is not an accepted 
theory. Results gives high values of {Q = O-Qto'l,) ([13])- ^ 

(e) Primordial nucleosynthesis: 
The theory of primordial nucleosynthesis, proposed in 1946 
by Gamov, assumes that the light elements till Lzy are gen- 
erated after big bang and that heavier elements originate 
from nuclear reactions inside stars. The values obtained 
for the abundances depends on some parameters like: rj, 
the value of the ratio baryons-photons, nowadays; 7V^, the 
number of neutrinos species; Tcmbr, the temperature of 
Cosmic Microwave Background Radiation. The theory of 
primordial nucleosynthesis permits to give limits to fJf, 
(b stands for baryons). With a ratio baryons-photons 
3 * 10-1° < 77 < 5 * 10^1", and a value of iV^, < 4 for 
the neutrinos species, Tcmrr — 2.736 ± O.OIK, is found 
0.011 < Slfe < 0.12 ([H]). 

In the following, we summarize the results of some more 
recent results. 

[l6j used several methods to obtain the value of O. Ac- 
cording to their classification, we divide the methods into 
the following four classes: 

• Global measures. Based on properties of space-time 
that constrain combinations of fi^ and the other cos- 
mological parameters (A, ffp, to). 

• Virialized Systems. Methods based on nonlinear 
dynamics within galaxies and clusters on comoving 
scales 1 - lOh-^Mpc. 

• Large-scale structure. Measurements based on mildly- 

nonlinear gravitational dynamics of fluctuations on 
scales 10 — 10Qh~^A'Ipc of supcrclustcrs and voids, in 
particular cosmic flows. 



• Growth rate of fluctuations. Comparisons of present 
day structure with fluctuations at the last scattering 
of the cosmic microwave background (CMB) or with 
high redshift objects of the young universe. 

The methods and current estimates are summarized in 
Table 3. The estimates based on virialized objects typi- 
cally yield low values of ^l,n ^ 0.2 — 0.3. The global mea- 
sures, large-scale structure and cosmic flows typically indi- 
cate higher values of ^Im ~ 0.4 — 1. 

Bahcall et al. ([l3l)) showed that the evolution of the 
number density of rich clusters of galaxies breaks the degen- 
eracy between J7 (the mass density ratio of the universe) and 
as (the normalization of the power spectrum), crgf2°'^ ~ 0.5, 
that follows from the observed present-day abundance of 
rich clusters. The evolution of high-mass (Coma-like) clus- 
ters is strong in f2 = 1, low-us models (such as the standard 
biased CDM model with erg ~ 0.5), where the number den- 
sity of clusters decreases by a factor of ~ 10'^ from z = 
to z ~ 0.5; the same clusters show only mild evolution in 
low-fi, high-fTg models, where the decrease is a factor of 
~ 10. This diagnostic provides a most powerful constraint 
on n. Using observations of clusters to z ~ 0.5 — 1, the 
authors found only mild evolution in the observed cluster 
abundance, and = 0.3 ± 0.1 and erg = 0.85 ± 0.15 (for 
A = models; for n + A = 1 models, = 0.34 ± 0.13). 

ferreira et al. ([lB|)j proposed an alternative method to 
estimate Vi2 directly from peculiar velocity samples, which 
contain redshift-independent distances as well as galaxy red- 
shifts. In contrast to other dynamical measures which de- 
termine (3 = fi^'^CTg, this method can provide an estimate of 
Jqj, range of cg where ft is the cosmological density 
parameter, while erg is the standard normalization for the 
power spectrum of density fluctuations. 

Melchiorri ( [l^ ) , used the angular power spectrum of the 
Cosmic Microwave Background, measured during the North 
American test flight of the BOOMERANG experiment, to 
constrain the geometry of the universe. Within the class of 
Cold Dark Matter models, they find that the overall frac- 
tional energy density of the universe, fl, is constrained to 
be 0.85 < O < 1.25 at the 68% confidence level. 

Branchini ([1^), compared the density and velocity fields 
as extracted from the Abell/ACO clusters to the corre- 
sponding fields recovered by the POTENT method from 
the Mark HI peculiar velocities of galaxies. Quantitative 
comparisons within a volume containing ^12 independent 
samples yield /3c = ^"-^/bc = 0.22 ± 0.08, where be is the 
cluster biasing parameter at 15h~^Mpc. If be ''^ 4.5, as in- 
dicated by the cluster correlation function, their result is 
consistent with 11 ~ 1. 

(f) Inflation: 

It is widely supposed that the very early universe experi- 
enced an era of inflation (see [2l[, [i3|). By 'inflation' 
one means that the scale factor has positive acceleration, 
a > 0, corresponding to repulsive gravity and 3p < —p. 
During inflation aH = d is increasing, so that comoving 
scales are leaving the horizon (Hubble distance) rather than 
entering it, and it is supposed that at the beginning of in- 
flation the observable universe was well within the horizon. 

The inflationary hypothesis is attractive because it holds 
out the possibility of calculating cosmological quantities, 
given the Lagrangian describing the fundamental interac- 
tions. The Standard Model, describing the interactions up 
to energies of order ITeV, is not viable in this context be- 
cause it does not permit inflation, but this should not be re- 
garded as a serious setback because it is universally agreed 
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mology. The nature of the required extension is not yet 
known, though it is conceivable that it could become known 
in the foreseeable future. But even without a specific model 
of the interactions (ie., a specific Lagrangian), the inflation- 
ary hypothesis can still offer guidance about what to expect 
in cosmology. More dramatically, one can turn around the 
theory-to-observation sequence, to rule out otherwise rea- 
sonable models. The importance of inflation is connected 
to: 

a) the origin of density perturbations, which could origi- 
nate during inflation as quantum fluctuations, which be- 
come classical as they leave the horizon and remain so on 
re-entry. The original quantum fluctuations are of exactly 
the same type as those of the electromagnetic field, which 
give rise to the experimentally observed Casimir effect. 

b) One of the most dramatic and simple effects is that there 
is no fine-tuning of the initial value of the density parame- 
ter ft = 8Trp/3m^iH^ . From the Friedmann equation, fl is 
given by 



(12) 



Its present value Qq is certainly within an order of mag- 
nitude of 1, and in the absence of an inflationary era fl 
becomes ever smaller as one goes back in time, implying an 
initial fine tuning. In contrast, if there is an inflationary era 
beginning when the observable universe is within the hori- 
zon, Eq. (fT2|) implies that flo will be of order 1, provided 
only that the same is true of ft at the beginning of infla- 
tion. A value of flo extremely close to 1 is the most natural, 
though it is not mandatory.'^ 

c) Another effect of inflation is that it can eliminate particles 
and topological defects which would otherwise be present. 
Anything produced before inflation is diluted away, and af- 
ter inflation there is a maximum temperature (the 'reheat' 
temperature) which is not high enough to produce all the 
particles and defects that might otherwise be present. As 
we shall remark later, this mechanism can remove desirable, 
as well as undesirable, objects. 

d) The most dramatic effect of inflation is that it may offer 
a way of understanding the homogeneity and isotropy of the 
universe, or at any rate of significant regions of it. We have 
nothing to say about this complex issue in its full generality, 
but a more modest version of it is our central concern. In 
this version, one begins the discussion at some early stage 
of inflation, when the universe is supposed already to be 
approximately homogeneous and isotropic. One then argues 
that in that case, scales far inside the horizon must be ab- 
solutely homogeneous and isotropic, except for the effect of 
vacuum fluctuations in the fields. Finally, one shows that 
after they leave the horizon, such fluctuations can become 
the classical perturbations that one deals with in cosmolog- 
ical perturbation theory. This pos sibility was first pointed 
out for gravitational waves by f23j and for density pertur- 
bations by several people ([24,]; [25i]; HI]). As we shall go to 
some trouble to demonstrate, the vacuum fluctuations can 
be evaluated unambiguously once an inflationary model is 
specified. 

(g) Scalar field inflation: 

Two mechanisms for inflation have been proposed. The 



^ An argument has been given for fig very close to 1 on the basis of 
effects on the cmb anisotropy from regions far outside the observable 



simplest one ([21[) invokes a scalar field, termed the infla- 
ton field. An alternative ([11]) is to invoke a modification 
of Einstein gravity, and combinations of the two mecha- 
nisms have also been proposed. During inflation however, 
the proposed modifications of gravity can be abolished by 
redefining the spacetime metric tensor, so that one recovers 
the scalar field case. We focus on it for the moment, but 
modified gravity models will be included later in our survey 
of specific models. 

In comoving coordinates a homogeneous scalar field 4> 
with minimal coupling to gravity has the equation of motion 



(j) + 3H(l3 + V'{(j)) = 
Its energy density and pressure are 

P = v + \^^ 



(13) 

(14) 
(15) 



If such a field dominates p and p, the inflationary condition 
3p < p is achieved provided that the field rolls sufficiently 
slowly. 



< V 



(16) 



Practically all of the usually considered models of inflation 
satisfy three conditions. First, the motion of the field is 
overdamped, so that the 'force' V' balances the 'friction 

term' 3H(j), 



Second, 



3H 



Wtt V V 



(17) 



(18) 



which means that the inflationary requirement (j)^ < V is 
well satisfied and 



1 87T 

3^ 



-V 



(19) 



These two conditions imply that H is slowly varying, and 
that the scale factor increases more or less exponentially, 



The third condition that is usually satisfied is 



M « 1 



where 



V 



V" 



(20) 



(21) 



(22) 



It can be 'derived' from the other two by differentiating 
the approximation Eq. ([TT]) for </) and noting that consis- 
tency with the exact expression Eq. (fT5)) requires ^ <C 1^' 
is satisfied. However there is no logical necessity for the 
derivative of an approximation to be itself a valid approx- 
imation, so this third condition is logically independent of 
the others. Conditions involving higher derivatives of V 
could be 'derived' by further differentiation, with the same 
caveat, but the two that we have given, involving only the 
first and second derivatives, are the ones needed to obtain 
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in which they are satisfied and we are adopting that nomen- 
clature here. Practically all of the usually considered models 
of inflation satisfy the slow-roll conditions more or less well. 

It should be noted that the first slow-roll condition is on 
a quite different footing from the other two, being a state- 
ment about the solution of the field equation as opposed to 
a statement about the potential that defines this equation. 
What we are saying is that in the usually considered models 
one can show that the first condition is an attractor solu- 
tion, in a regime typically characterized by the other two 
conditions, and that moreover reasonable initial conditions 
on will ensure that this solution is achieved well before 
the observable universe leaves the horizon. It is important 
to remember that there are strong observational limits for 
the parameters previously introduced (e.g. e, rj). For ex- 
ample [l^l studied the possible contribution of a stochastic 
gravitational wave background to the anisotropy of the cos- 
mic microwave background in cold and mixed dark matter 
(CDM and MDM) models. This contribution was tested 
against detections of CMB anisotropy at large and inter- 
mediate angular scales. The best fit parameters (i.e. those 
which maximize the likelihood) are (with 95% confidence) 
ns 



1.23l^;i^ and 



29e 



2 4+3 '* 



where 



fins) 



r(3-ns)r( — - — ) 

r2(i:i:^)r(^^^ 

^ 2 ^ ^ 2 



(23) 



(24) 



The previous constraint fixes the value of e as well that of i] 



21] = Us ~ 1 + 2e 



(25) 



They find that by including the possibility of such back- 
ground in CMB data analysis it can drastically alter the 
conclusion on the remaining cosmological parameters. More 
stringent constraints on some of the previous parameters are 
given in section 1.12. 

(h) Conclusions: 
We have seen the possible values of il using different meth- 
ods. We have to add that Cosmologists are "attracted" by 
a value of ilo = 1. This value of is requested by infla- 
tionary theory. The previous data lead us to the following 
hypotheses: 

i) r^o < 0.12; in this case one can suppose that the uni- 
verse is fundamentally made of baryonic matter (black holes; 
Jupiters; white dwarfs). 

ii) fio > 0.12; in this case in order to have a fiat universe, it 
is necessary a non- baryonic component, fib = 1 is excluded 
by several reasons (see [1^ , [l3| • The remaining possibilities 
are: 

1) existence of a smooth component with ft = 0.8. 

The test of a smooth component can be done with kinematic 
methods. 

2) Existence of a cosmological term, absolutely smooth to 
whom correspond an energy density p^ac = 

3) existence of non-baryonic matter: the universe is fun- 
damentally done of particles (neutrinos, WIMPS (Weakly 
Interacting Massive Particles)). 

4) A combination of 2) and 3). 

Before going one, I want to recall that some authors ([29*]; 
(30| : [3l| ) have assumed that we have a scant knowledge 
of physical laws. Sanders assumes that the gravitational 



Milgrom assumes that the Newton law of gravitation is not 
valid when the gradient of the potential is small. In this 
case, the problem of the dynamics of clusters of galaxies is 
solved without introducing Dark Matter. In any case, the 
quoted assumptions have no general theory that can justify 
them. 

1.2 Dark matter in particles 

We know that if = 1 , dark matter cannot be constituted 
exclusive of baryonic matter. The most widespread hypoth- 
esis is that dark matter is in form of particles. Several candi- 
dates exist: neutrinos, axions, neutralinos, photinos, grav- 
itinos, etc. Interesting particles are usually grouped into 
three families: 

HDM (Hot Dark Matter), CDM (Cold Dark Matter) and 
WDM (Warm Dark Matter). In order to understand this 
classification it is necessary to go back to the early phases 
of universe evolution. The history of the universe is charac- 
terized by long phases of local thermodynamic equilibrium 
(LTE) and by "deviation" by it: nucleosynthesys, bario- 
genesys, decoupling of species, etc. In the early universe, 
were present the particles that we know today and other 
particles predicted theoretically, but that have not been 
observed. Massive particles preserved the thermodynamic 
equilibrium concentration until the rate, F, of reactions and 
interactions that produced that concentration was larger 
than the expansion rate of the universe, H. When the con- 
dition r > H was no longer satisfied the reactions stopped 
and the abundance of the considered species remained con- 
stant at the value it had at time of freez-out, Tf, time at 
which T — H. If we indicate with Y = - the number of 
particles per unit comoving volume and we remember that 
n is the number density of species and s the entropy density, 
we obtain a contribution of the species to the actual density 
of the universe as Qh^ = 0.28Y{Tf){^) (^). At time T/, 
particles could be relativistic or non-relativistic. Relativis- 
tic particles are today indicated with the term hot cosmic 
relics, HDM, while non-relativistic particles are named cold 
cosmic relics, CDM. There is an intermediate case, that of 
warm relics, WDM. 

An example of HDM are massive neutrinos. Possible masses 
for these neutrinos are: 



25ew < nil, < lOOew 



(26) 



li) and 



4.9— 1.3Gei; for Maiorana's neutrinos 

1.3 — 4.2Gew for Dirac's neutrinos ^ ' 



There are confirmed experimental evidence of the existence 
of massive neutrinos. In 1980, ^33] announced the detection 
of an electronic antineutrino with mass 30 eV, by means 
of the shape of the electron energy spectrum in the f3 
decay of tritium ('si']). Experiments (Super-Kamiokande, 
SNOW) have obtained some evidence of non-zero mass 
from neutrino oscillations. This yields a difference of square 
masses of order 10"'^ eV, and a mass of 0.05 eV (in the 
simplest case) (see [35|). 

Among typical examples of CDM we have WIMPS and in 
particular axions and neutralinos (SUSY particle). This 
particle was postulated in order to solve the strong CP 
problem in nuclear physics. This problem arises from the 
fact that some interactions violate the parity, P, time inver- 
sion, T, and CP. If these are not eliminated, they give rise 
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limits (pS'l)- The solution to the problem was proposed by 
Peccei-Quinn in 1977 ([s^) in terms of a spontaneous sym- 
metry breaking scheme. To this symmetry breaking should 
be associated a Nambu-Goldstone boson: the axion. The 
axion mass ranges between 10~^^ev-lGev. In cosmology 
there are two ranges of interest: 10~^eti < nia < 10~^ei; ; 
Sev < rUa < Sew. Axion production in the quoted range can 
originate due to a series of astrophysical processes f[l3l|) 
and several are the ways these particles can be detected. 
Nevertheless the effort of researchers expecially in USA, 
Japan and Italy, axions remain hypothetical particles. 
They are in any case the most important CDM candidates. 

In the following, I am going to speak about the basic 
ideas of structure formation. I shall write about density 
perturbations, their spectrum and evolution, about correla- 
tion functions and their time evolution, etc. 

1.3 Origin of structures 

Observing our universe, we notice a clear evidence of in- 
homogeneity when we consider small scales (Mpcs). In clus- 
ters density reaches values of 10'^ times larger than the av- 
erage density, and in galaxies it has values 10^ larger than 
the average density ([i3(). If we consider scales larger than 
10^ Mpcs universe appears isoptric as it is observed in the 
radio-galaxies counts, in CMBR, in the X background (jTljl). 
The isotropy at the decoupling time, tdec, at which matter 
and radiation decoupled, universe was very homogeneous, 
as showed by the simple relation: 



dp ST 
— — const—— 
P T 



(28) 



([3) ^- The difference between the actual universe and 
that at decoupling is evident. The transformation between 
a highly homogeneous universe, at early times, to an highly 
local non homogeneous one, can be explained supposing that 
at tdec were present small inhomogeneities which grow up 
because of the gravitational instability mechanism (p37jV 
Events leading to structure formation can be enumerated 
as follows: 

(a) Origin of quantum fluctuations at Planck epoch. 

(b) Fluctuations enter the horizon and they grow linearly 
till recombination. 

(c) Perturbations grow up in a different way for HDM and 
CDM in the post-recombination phase, till they reach the 
non-linear phase. 

(d) Collapse and structure formation. 

Before tdec inhomogeneities in baryonic components could 
not grow because photons and baryons were strictly cou- 
pled. This problem was not present for the CDM compo- 
nent. Then CDM perturbations started to grow up before 
those in the baryonic component when universe was mat- 
ter dominated. The epoch teq ~ 4.4 * 10^"(rio^^)~^sec, at 
which matter and radiation density are almost equal, can 
be considered as the epoch at which structures started to 
form. The study of structure formation is fundamentally an 
initial value problem. Data necessary for starting this study 
are: 

1) Value of r^Q. In CDM models the value chosen for this 
parameter is 1, in conformity with inflationary theory pre- 
dictions. 

2) The values of Q,i for the different components in the uni- 
verse. For example in the case of baryons, nucleosynthesis 



gives us the limit 0.014 < fib < 0.15 while ^wimps ~ 0.9. 
3) The perturbation spectrum and the nature of pertur- 
bations (adiabatic or isocurvature). The spectrum gener- 
ally used is that of Harrison-Zeldovich: P{k) = Afc" with 
n — 1. The perturbation more used are adiabatic or curva- 
ture. This choice is dictated from the comparison between 
theory and observations of CMBR anisotropy. 

H The spectrum of density perturbation 

In order to study the distribution of matter density in 
the universe it is generally assumed that this distribution 
is given by the superposition of plane waves independently 
evolving, at least until they are in the linear regime (this 
means till the overdensity 5 — << 1). Let we divide 
universe in cells of volume Vu and let we impose periodic 
conditions on the surfaces. If we indicate with p the average 
density in the volume and with p[r) the density in r, it is 
possible to define the density contrast as: 



5{v) 



p{t^) - P 
P 



This quantity can be developed in Fourier series: 
S{r) — ^ S\^exp{i]i.r) — ^ S\^exp{~ik.r) 



(29) 



(30) 



([Tj]), where 



27m a- 



(and similar conditions for the other 



components) and for the periodicity condition S{x,y,L) = 
5{x, y, 0) (and similar conditions for the other components). 
Fourier coefficients 5k are complex quantities given by: 



5k = — / (5(r)ea;p(— ikr)(ir 



(31) 



For mass conservation in Vu we have also (5k=o = while for 
reality of i5(r), 5^ = 5_k- If we consider n volumes, Vu, we 
have the problem of determining the distribution of Fourier 
coefficients 5k and that of \5\. We know that the coefficients 
are complex quantities and then (5k = \Sk\ exp{i6k) . If we 
suppose that phases are random^ in the limit — > cx) it is 
possible to show that we get \S\ = l^kT- The Central 
limit theorem leads us to conclude that the distribution for 
6 is Gaussian: 



P{S) (X exp{ — 



(32) 



(|28|). The quantity a that is present in Eq. ((32)) is the 
variance of the density field and is defined as: 

a'=<S'>=J2<\Sk\'>^^Y.^l (33) 

^u 



This quantity characterizes the amplitude of the inhomo- 
geneity of the density field. If T4i — > oo, we obtain the more 
usual relation: 



(2^) 



k^dk 



(34) 



The term P{k) =< \5\^ > is called "Spectrum of perturba- 
tions" . It is function only of k because the ensemble average 
in an isotropic universe depends only on r. A choice often 
made for the primordial spectrum is P{k) = Ak^ which 
in the case n — \ gives the scale invariant spectrum of 
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It can be defined as the joint probabihty of finding an over- 
density 6 in two distinct points of space: 



C(r,t) ^(r,t)(5(r + x,i) > 



(35) 



(p3q). where averages are averages on an ensemble obtained 
from several realizations of universe. Correlation function 
can be expressed as the joint probability of finding a galaxy 
in a volume SVi and another in a volume SV2 separated by 
a distance ri2: 



(36) 



where ny is the average number of galaxies per unit volume. 
The concept of correlation function, given in this terms, can 
be enlarged to the case of three or more points. 
Correlation functions have a fundamental role in the study 
of clustering of matter. If we want to use this function for 
a complete description of clustering, one needs to know the 
correlation functions of order larger than two ([39|). By 
means of correlation functions it is possible to study the 
evolution of clustering. The correlation functions are, in 
fact, connected one another by means of an infinite system 
of equations obtained from moments of Boltzmann equa- 
tion which constitutes the BBGKY (Bogolyubov-Green- 
Kirkwood-Yvon) hierarchy ([40|). This hierarchy can be 
transformed into a closed system of equation using closure 
conditions. Solving the system one gets information on cor- 
relation functions. 

In order to show the relation between perturbation spec- 
trum and two-points correlation function, we introduce in 
Eq. ([351), Eq. recalling that = S{-k) and taking the 
limit Vu 00, the average in the Eq. ()35|1 can be expressed 
in terms of the integral: 



e(r)- 



1 



(27r)3 



\S{k)\'^exp{-ikr)d^k 



(37) 



This result shows that the two-point correlation function 
is the Fourier transform of the spectrum. In an isotropic 
universe, it is |r| = r and then |k| — k and the spectrum can 
be obtained from an integral on |k| = k. Then correlation 
function may be written as: 



1 /" , 9 „/, V sin(kr) ,, 
2^ / ^^^(^)^* 



(38) 



During the evolution of the universe and after perturbations 
enter the horizon, the spectrum is subject to modulations 
because of physical processes characteristic of the model it- 
self (Silk damping ([4l|) for acoUisional components, free 
streaming for coUisional particles, etc.). These effects are 
taken into account by means of the transfer function T(k; t) 
which connects the primordial spectrum P{k;tp) at time tp 
to the final time tj: 



P{k-tf) 



Htf) 



h{tp) 



T\k-tf)P{k-tp) 



(39) 



where b(t) is the law of grow of perturbations, in the linear 
regime. In the case of CDM models the transfer function is: 

-1 

T{k) = |l + [afc -f (6fc)i-5 + [ckfY] " (40) 

(Hi), where a = ^AiSlh^Y^Mpc; b = 3.0{nh'^)-Hlpc; c = 
1.7{flh^)^^ Mpc; V — 1.13. It is interesting to note that Eq. 
([32|) is valid only if cr << 1, since \5\ < 1. This implies than 



certain point modes are no longer independent and start 
to couple giving rise to non-linear effects that change the 
spectrum and correlation function (|43]). There are also 
some theories (e.g., cosmic strings (j44i])) in which even in 
the linear regime perturbations are not Gaussian. 

1.5 Curvature and isocurvature perturbations 

The study of the evolution of density perturbations can 
be divided into two phases: 

1) perturbations are outside the horizon, in other terms 
they have a scale A larger than Hubble radius rn — ct or 
A > H-i; 

2) perturbations are inside the horizon, A < H^^. 

In studying the first case it is necessary to use general rela- 
tivity and one can demonstrate that two different kinds of 
fiuctuations exist: curvature or adiabatic (motivated by the 
simplest models of inflation) and isocurvature or isothermal. 
Curvature fluctuations are characterized by a fluctuation in 
energy density or in space curvature. For them, we can 
write: 



5S_ 



3 5pr 

4 pr 



6T 

> 



Spr, 
Prn 







(41) 



f|28l|). where with m has been indicated the matter com- 
ponent, with r radiation, with S entropy and with T tem- 
perature. Last equation explains why these fluctuations are 
named adiabatic, since for them the entropy variation is 
zero. 

Isocurvature or isothermal perturbations are not character- 
ized by fluctuations in the curvature of the metric, but they 
are fluctuations in the local equation of state of universe and 
in agreement with the name it results ST = 0. Until fluctua- 
tions of isocurvature are not inside the horizon, the causality 
principle does not permit a redistribution of energy density. 
This is possible only when perturbations enter the Horizon 
and isocurvature perturbations can be converted into per- 
turbations in the energy density. As a consequence, the 
distinction between those two kinds of perturbations is no 
longer meaningful after them enter the horizon ([4^; 0|)- 
The origin of curvature perturbations may be explained in- 
side inflationary models or assuming that they are initially 
present as perturbations of the metric. Isocurvature fluc- 
tuations may be always produced in inflationary scenarios 
from fluctuations in the density number of barions or axions 

(Hi, [11). 

1.6 Perturbations evolution 

Density perturbations in the components of the universe 
evolve with time. In order to get the evolution equations for 
S in Newtonian regime, it is possible to use several models. 
In our model, we assume that gravitation dominates on the 
other interactions and that particles (representing galaxies, 
etc.) move coUisionless in the potential ^ of a smooth den- 
sity function ([38!]). The distribution function of particles 
for position and momentum is given by: 



dN = f{x,p,t)d^xd^p 



and density: 
p(x,t) = 



,-3 



J d^pf{x,p,t)=pt, ll + 5(x,i)] 



(42) 



(43) 



where m is the mass of a particle and pb the background 
density. Applying Liouville theorem to the probability den- 
sity on a limited region of phase-space of the system we have 
that f verifies the equation: 
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The distribution function f that appears in the previous 
equations cannot be obtained from observations. It is possi- 
ble to measure moments off (density, average velocity, etc.). 
We want now to obtain the evolution equations for 5. For 
this reason, we start integrating Eq. (|44|) on p and after 
using Eq. (gS]), we get: 



If we define velocity as: 

^ J ma ^ 

Jfd^p 

and introduce it in Eq. (I45p we get: 

Pb^ + - V = 
ot a 



(45) 



(46) 



(47) 



We can now multiply Eq. (l44l) for p and integrate it on the 
momentum: 

^ / Vo.fd'p+^dp f p^ppfd?p+a:'p{^,t)4>,^^0 (48) 
ot / ma^ / 



this last in Eq. leaves us with: 



adS 
adt 



:V[(1 + '^) V 



1 



PbTTia' 



and finally using 



dad/3 I PaP0(f><fp 

(49) 



< VaVf3 > = 



/ fPaPpd^P 



(50) 



ma^ J f(Pp 

the equation for the evolution of overdensity becomes: 

^ + ^) ^ + [(1 + 5)< v-v^ >] 

_ (51) 
f[38|V The term < VaVp > is the tensor of anisotropy of 
peculiar velocity. This is present in the gradient and then it 
behaves as a pressure force. If we consider an isolated and 
spherical perturbation, it is possible to assume that initial 
asymmetries does not grow up and so we can suppose, in 
this hypothesis that < VaV/3 >= 0. In this case and with 
the linearity assumption 6 << 1 we have: 



at'' a ot 



(52) 



This equation in an Einstein-de Sitter universe (fi = 1, A 
0) has the solutions: 



5+=A+(x)t = 



(5_(x,i) = A_(x)i-i 



(53) 



The perturbation is then done of two parts: a growing one, 
becoming more and more important with time, and a de- 
caying one becoming negligible with increasing time, with 
respect to the growing one. 

In the case of open models with no cosmological constant: 

< 1, A = 0, we can write: 



(54) 



0(77) 



no 



2(1 -r!o) 



2i7o(l-f^o)='/2 



(coshyy — 1) (55) 
(sinhyy — rj). 



In the case of flat models with positive cosmological con- 
stant: n < 1, A 0, Q + A/3Hq = 1, we can write: 



^^2^nGp{i + {n,'-i) 



a{t) = {n^- 



l)"'/%inh2/3 




(56) 



(57) 



Before concluding this section, we want to find an ex- 
pression for the velocity field in the linear regime. Using 
the equation of motion p ~ ma^i, ^ — —m y (j) and the 
proper velocity of a particle, v — ax, verify the equation: 



^ + = = Gp,a I S.8i..t)^^ (58) 

at a a J |x — x | 

Supposing that v is a similar solution for the density, v = 
V+(x,t)tP, we get: 



47r dx° 



(59) 



(|38l|). This solution is valid just as that for 5 in the linear 
regime. At time t — to this regime is valid on scales larger 
than Sh^'^Mpc. 

1. 7 Non-linear phase 

Linear evolution is valid only if i5 << 1 or similarly, if the 
mass variance, cr, is much less than unity. When this condi- 
tion is no longer verified (e.g., if we consider scales smaller 
than 8/1^^ Mpc), it is necessary to develope a non-linear 
theory. In regions smaller than 8/1^^ Mpc galaxies are not a 
Poisson distribution but they tend to cluster. If one wants 
to study the properties of galactic structures or clusters of 
galaxies, it is necessary to introduce a non-linear theory of 
clustering. A theory of this last item is too complicated to 
be developed in a purely theoretical fashion. The problem 
can be faced assuming certain approximations that simpli- 
fies it ([13]) or as often it is done, by using N-Body simu- 
lations of the interesting system. The approximations are 
often used to furnish the initial data to simulations. In the 
simulations, a large number of particles are randomly dis- 
tributed in a sphere, in the points of a cubic grid, in order 
to eliminate small scale noise. The initial spectrum is ob- 
tained perturbing the initial positions by means of a super- 
position of plane waves having random distributed phases 
and wave vector ([48^1). Obviously, the universe is considered 
in expansion (or comoving coordinates are used), and then 
the equation of motion of particles are numerically solved. 
For what concerns the analytical approximations one of the 
most used is that of 47]. This gives a solution to the prob- 
lem of the grow of perturbations in an universe with p = 
not only in the linear regime but even in the mildly non- 
linear regime. In this approximation, one supposes to have 
particles with initial position given in Lagrangian coordi- 
nates q. The positions of particles, at a given time t, are 
given by: 
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where x indicates the Eulerian coordinates, p(q) describes 
the initial density fluctuations and b{t) describes their grow 
in the hnear phase and it satisfies the equation: 



— ^ .dh da . „ , 
-'"^ dt dt ' 



(61) 



The equation of motion of particles, according to the quoted 
approximation, is given by: 



V = aq + 6p(q) 
The peculiar velocity of particles is given by: 

dx db 

while the density of the perturbed system is given by: 



p(q,0 





dqj 




p 


dxk 


= P 



6(t)- 



(62) 



(63) 



(64) 



Developing the Jacobian present in Eq. (|64p at first order 
in 5(t)p(q), one obtains: 



Sp 
P 



-b{t) VqP(q) 



(65) 



This equation can be re-written, separating the space and 
time dependence, as in the equation for u, and writing: 



h{t) = ti 
in the form: 



P(q) ^^i--^Ai,exp{ik.q) 



5p 



— = ^ Ak<3ea;p(ikq) 



P 



(66) 



(67) 



([28|), that leads us back to the linear theory. In other 
words, Ze'ldovich approximation is able to reproduce the 
linear theory, and is also able to give a good approximation 
in regions with ^ >> 1. Using the expression for p{q), the 



Jacobian in Eq. (j64p is a real matrix and symmetric that 
can be diagonalized. With this p{q) the perturbed density 
can be written as: 

^^"^'^^ = (1 - 6(t)Ai(q))(l - h{t)\^{q)){\ - h{t)M<l)) 

(68) 

where Ai, A2, A3 are the three eigenvalues of the Jacobian, 
describing the expansion and contraction of mass along the 
principal axes. From the structure of the last equation, we 
notice that in regions of high density Eq. (|68p becomes 
infinite and the structure of collapse in a pancake, in a 
filamentary structure or in a node, according to values of 
eigenvalues. Some N-body simulations ([i^) tried to ver- 
ify the prediction of Ze'ldovich approximation, using initial 
conditions generated using a spectrum with a cut-off at low 
frequencies. The results showed a good agreement between 
theory and simulations, for the initial phases of the evolu- 
tion {a{t) — 3.6). Going on, the approximation is no more 
valid starting from the time of shell-crossing. After shell- 
crossing, particles does not oscillate any longer around the 
structure but they pass through it making it vanish. This 
problem has been partly solved supposing that particles, 
before reaching the singularity they sticks the one on the 
other, due to a dissipative term that simulates gravity and 
then collects on the forming structure. This model is known 
as "adesion-model" (jSOj). 

Summarizing, Zel'dovich approximation gives a description 
of the transition between linear and non-linear phase. It 



1.8 Quasi-linear regime 

We have seen in the previous section that in the case 
of regions of dimension smaller than 8h~^AIpc, the linear 
theory is no more a good approximation and a new theory 
is needed or N-body simulations. Non-linear theory is able 
to calculate quantities as the formation redshift of a given 
class of objects as galaxies and clusters, the number of 
bound objects having masses larger than a given one, the 
average virial velocity and the correlation function. It is 
possible to get an estimate of the given quantities as that 
of other not cited, using an intermediate theory between 
the linear and non-linear theory: the quasi-linear theory. 
This last is obtained adding to the linear theory a model 
of gravitational collapse, just as the spherical collapse 
model. Important results that the theory gives is the 
bottom- up formation of structures (in the CDM model). 
Other important results are obtained if we identify density 
peaks in linear regime with sites of structure formation. 
Two important papers in the development of this theory 
are [5l| and that of [s^l- This last paper is an application 
of the ideas of the quasi-linear theory to the CDM model. 
The principles of this approach are the following: 



• Regions of mass larger than M that collapsed can 
be identified with regions where the density contrast 
evolved according to linear regime, S{M,x), has a 
value larger than a threshold, 6c. 

• After collapse regions does not fragment. 

The major drawbacks of the theory, as described in [s^ are 
fundamentally the fact that the estimates that can be ob- 
tained by means of this theory depends on the threshold 
Sc, on the ratio between the filtering mass and that of ob- 
jects and from other parameters. Nevertheless, this theory 
has helped cosmologists in obtaining estimate of important 
quantities as those previously quoted, and at same time give 
evidences that leads to exclude very low values for spectrum 
normalization. 

1.9 Spherical Collapse 

Spherical symmetry is one of the few cases in which grav- 
itational collapse can be solved exactly ([11]; [1^). In fact, 
as a consequence of Birkhoff's theorem, a spherical pertur- 
bation evolves as a FRW Universe with density equal to the 
mean density inside the perturbation. 

The simplest spherical perturbation is the top-hat one, 
i.e. a constant overdensity S inside a sphere of radius R; 
to avoid a feedback reaction on the background model, the 
overdensity has to be surrounded by a spherical underdense 
shell, such to make the total perturbation vanish. The evo- 
lution of the radius of the perturbation is then given by a 
Friedmann equation. 

The evolution of a spherical perturbation depends only on 
its initial overdensity. In an Einstein-de Sitter background, 
any spherical overdensity reaches a singularity (collapse) at 
a final time: 



3^ /5,, A"^/' 



(69) 



By that time its linear density contrast reaches the value: 



0/0 \ 



3/2 
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In an open Universe not any overdensity is going to col- 
lapse: the initial density contrast has to be such that 
the total density inside the perturbation overcomes the 
critical density. This can be quantified (not exactly but 
very accurately) as follows: the growing mode saturates at 
b{t) = 5/2($1q ^ — 1), so that a perturbation ought to satisfy 
Si > 1.69 • 2{n^^ - l)/5 to be able collapse. 

Of course, collapse to a singularity is not what really hap- 
pens in reality. It is typically supposed that the structure 
reaches virial equilibrium at that time. In this case, argu- 
ments based on the virial theorem and on energy conserva- 
tion show that the structure reaches a radius equal to half 
its maximum expansion radius, and a density contrast of 
about 178. In the subsequent evolution the radius and the 
physical density of the virialized structure remains constant, 
and its density contrast grows with time, as the background 
density decays. Similarly, structures which collapse before 
are denser than the ones which collapse later. 

Spherical collapse is not a realistic description of the 
formation of real structures; however, it has been shown 
(see [M] for a rigorous proof or [sl], [s^) that high peaks 
(> 2a) follow spherical collapse, at least in the first phases 
of their evolution. However, a small systematic departure 
from spherical collapse can change the statistical properties 
of collapse times. 

Spherical collapse can describe the evolution of under- 
densities. A spherical underdensity is not able to collapse 
(unless the Universe is closed!), but behaves as an open 
Universe, always expanding unless its borders collide with 
neighboring regions. At variance with overdensities, under- 
densities tend to be more spherical as they evolve, so that 
the spherical model provides a very good approximation for 
their evolution. 

1.9.1 Improvements to the Spherical Collapse model 

Several years ago it was realized that the density field dis- 
tributions around the density peaks, which eventually will 
give birth to galaxies and clusters, depart from spherical 
symmetry and from the average density profile, producing 
important consequences on collapse dynamics and forma- 
tion of protostructures ([13; [11; [Hi]; [H^). A 
fundamental role in this context is played by the joint ac- 
tion of tidal torques (coupling shells of matter which are 
accreted around a density peak and neighbori ng p rotostruc- 
tures (dll; H^)), and by dynamical friction ([1^; [13], [HI, 

(H). 

According to the previrialization conjecture ([H^, [6^). 
initial asphericities and tidal interactions between neighbor- 
ing density fluctuations induce significant non-radial mo- 
tions which oppose the collapse. This means that virialized 
clumps form later, with respect to the predictions of the 
linear perturbation theory or the spherical collapse model 
(hereafter SM), and that the initial density contrast, needed 
to obtain a given final density contrast, must be larger than 
that for an isolated spherical fluctuation. This kind of con- 
clusion was supported by [66j , [63| , , [1^ and [T^I . 
In particular [6a] and [63| pointed out that non-radial mo- 
tions would slow the rate of growth of the density con- 
trast by lowering the peculiar velocity and suppress collapse 
once the system detaches from general expansion. [6^ gave 
examples of the growth of non-radial motions in N-body 
simulations. Arguments based on a numerical least-action 
method lead 64] to the conclusion that irregularities in the 
mass distribution, together with external tides, induce non- 
radial motions that slow down the collapse. [t^I used N- 



the slope of the initial power spectrum is n > —1, non-linear 
tidal interactions slow down the growth of density fluctua- 
tions and the magnitude of the effect increases when n is 
increased. 

Opposite conclusions were obtained by [lH, [ll], [t^I, 
In particular (Tlj . using the quasi- linear (QL) approxima- 
tion ([131) showed that the shear affects the dynamics of 
collapsing objects and it leads to infall velocities that are 
larger than in the case of non-shearing ones. Bertschinger 
& Jain {.73\) put this result in theorem form, according to 
which spherical perturbations are the slowest in collapsing. 
Bartelmann et al. (1993) argued that the collapse does not 
start from a comoving motion of the perturbation, but that 
the continuity equation requires an initial velocity perturba- 
tion directly related to the density perturbation. The effect 
is that collapse proceeds faster than in the case where the 
initial velocity perturbation is set equal to zero and the col- 
lapse timescale is shortened. The N-body simulations by 



72l | did not reproduce previrialization effect, but the reason 
is due to the fact that they assumed an n = — 1 spectrum, 
differently from the n = one used by [13] that reproduced 
the effect. If n < —1 the peculiar gravitational acceleration, 
g oc diverges at large R and the gravitational 

acceleration moves the fluid more or less uniformly, gener- 
ating bulk flows rather than shearing motions. Therefore, 
its collapse will be similar to that of an isolated spherical 
clump. If 71 > —1, the dominant sources of acceleration 
are local, small-scale inhomogeneities and tidal effects will 
tend to generate non-radial motions and resist gravitational 
collapse. In a more recent paper, (75| have proposed some 
analytic prescriptions to compute the collapse time along 
the second and the third principal axes of an ellipsoid, by 
means of the 'fuzzy' threshold approach. They pointed out 
that the formation of real virialized clumps must correspond 
to the third axis collapse and that the collapse along this 
axis is slowed down by the effect of the shear rather than 
be accelerated by it, in contrast to its effect on the first 
axis collapse. They concluded that spherical collapse is the 
fastest, in disagreement with Bertschinger & Jain's theo- 
rem. This result is in agreement with [62i] . The quoted 
controversy was addressed by [t^ who examined the evolu- 
tion of non-spherical inhomogeneities in a Einstein-dc Sitter 
universe, by numerically solving the equations of motion for 
the principal axes and the density of a dust ellipsoid. They 
showed that for lower values oi v {v = 2) the growth rate 
enhancement of the density contrast induced by the shear 
is counterbalanced by the effect of angular momentum ac- 
quisition. For > 3 the effect of angular momentum and 
shear reduces, and the evolution of perturbations tends to 
follow the behavior obtained in the SM. [tJI studied the 
role of shear fields on the evolution of density perturbations 
by using an analytical approximate solution for the equa- 
tions of motion of homogeneous ellipsoids embedded in a 
homogeneous background. The equations of motion of a 
homogeneous ellipsoid ([ll; [li] (hereafter WS)) were mod- 
ified in order to take account of the tidal field, as done in 
(soj and then were integrated analytically, similar to what 
was done in WS. The density contrast at turn-around and 
the collapse velocity were found to be reduced with respect 
to that found by means of the SM. The reduction increases 
with increasing strength of the external tidal field and with 
increasing initial asymmetry of the ellipsoids. 

The second physical effect with changes cluster collapse 
is dynamical friction. Former treatments of the dynami- 
cal friction effects on the structure of clusters of galaxies, 
considering only the component generated by the galactic 
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namical friction taking into account the effect of substruc- 
ture, showing that dynamical friction delays the collapse of 
low-v peaks inducing a bias of dynamical nature. Because 
of dynamical friction under-dense regions in clusters (the 
clusters outskirts) accrete less mass with respect to that 
accreted in absence of this dissipative effect and as a con- 
sequence over-dense regions are biased toward higher mass 
(AC). Dynamical friction and non-radial motions acts in a 
similar fashion: they delay the shell collapse consequently 
inducing a dynamical bias. Whenever efficient, these mech- 
anisms will generate a physical selection of those peaks in 
the initial density field that eventually will give rise to the 
observed cosmic structures. As a consequence of dynami- 
cal friction and tidal torques, one expects changes in the 
threshold of collapse, the mass function and the correlation 
function. 

In the next subsections, I shall study how the spherical 
collapse model is changed by the joint effect of dynamical 
friction, tidal torques and a non-zero cosmological constant. 

1.9.2 Dynamical friction and structure formation 

In a hierarchical structure formation model, the large 
scale cosmic environment can be represented as a collision- 
less medium made of a hierarchy of density fluctuations 
whose mass, M, is given by the mass function N{M,z), 
where z is the redshift. In these models matter is con- 
centrated in lumps, and the lumps into groups and so on. 
In such a material system, gravitational field can be de- 
composed into an average field, Fo(r), generated from the 
smoothed out distribution of mass, and a stochastic compo- 
nent, Fstoch{f), generated from the fluctuations in number 
of the field particles. The stochastic component of the grav- 
itational field is specified assigning a probability density, 
W{F), ([82]). In an infinite homogeneous unclustered sys- 
tem iy(F) is given by Holtsmark distribution ([l^l) while 
in inhomogeneous and clustered systems VF(F) is given by 
[S^ l and [U respectively. The stochastic force, F stock, in a 
self-gravitating system modifies the motion of particles as it 
is done by a frictional force. In fact a particle moving faster 
than its neighbors produces a defiection of their orbits in 
such a way that average density is greater in the direction 
opposite to that of traveling causing a slowing down in its 
motion. Following [s^] method, the frictional force which is 
experienced by a body of mass M (galaxy), moving through 
a homogeneous and isotropic distribution of lighter parti- 
cles of mass m (substructure), having a velocity distribution 
n{v) is given by: 



M- 



dt 



= -AiiG^M'^n{v)^ log Ap 



(71) 



where log A is the Coulomb logarithm, p the density of the 
field particles (substructure). 

A more general formula is that given by [ssl ] in the hypoth- 
esis that there are no correlations among random force and 
their derivatives: 

f W(F)F^T(F)d^F , , 

2 < > 

where 77 is the coefiicient of dynamical friction, T(i^) the 
average duration of a random force impulse, < > the 
characteristic speed of a field particle having a distance r ~ 
(QM_Yl'^ from a test particle (galaxy). This formula is more 
general than Eq. (|7T|) because the frictional force can be 
calculated also for inhomogeneous systems when W(F) is 
given. If the field particles are distributed homogeneously 
the dynamical friction force is given by: 



f|83l|). where rua and ria are respectively the average mass 
and density of the field particles. Using virial theorem we 
also have: 



1/3 



m n 



(74) 



sys 



where Mtot is the total mass of the system, Rsys its radius 
and N is the total number of field particles. The dynamical 
friction force can be written as follows: 



F — —rjv ~ — 



A.U[Gfnanac 



1/2 



N 



where N = ^Kysn^ 



(75) 

^ is the comoving 



and riac = x a' 
number density of peaks of substructure of field particles. 
This last equation supposes that the field particles gener- 
ating the stochastic field are virialized. This is justified by 
the previrialization hypothesis ( 65]). 

To calculate the dynamical evolution of the galactic com- 
ponent of the cluster it is necessary to calculate the num- 
ber and average mass of the field particles generating the 
stochastic field. 

The protocluster, before the ultimate collapse at z ~ 0.02, 
is made of substructure having masses ranging from 10^ — 
IO^Mq and from galaxies. I suppose that the stochastic 
gravitational field is generated from that portion of sub- 
structure having a central height v larger than a criti- 
cal threshold Vc. This latter quantity can be calculated 
(following AC) using the condition that the peak radius, 
^pk{i^ > Vc), is much less than the average peak separation 
TT-aiv > Vc)^^/'^ where iia is given by the formula of [52l| for 
the upcrossing points: 



nac{v > Vc) 



^ exp(i.g/2) 7 3 , 
(2^)2 ^rJ ^'^^ 

4V3 ... 2 2 



1 



572(1-572/9)1/2 



exp(-57'^.,7l8)] 



(76) 



where 7, i?* are parameters related to moments of the power 
spectrum ([52] Eq. 4.6A). The condition rpk{v > Vc) < 
Q.lriaiy > Vc)^^/^ ensures that the peaks of substructure 
are point like. Using the radius for a peak: 



rpk 



V2R^ 



-,1/3 



(l + i^ao)(73 + (0.9H)3/2 



(77) 



(AC), I obtain a value of Vc — 1.3 and then we have na{v > 
Vc) — 50.7Mpc~'^ (7 = 0.4, R^ = bOKpc) and nia is given 
by: 



na{v > Vc) 



mpkiv)Npk{v)diy = 10^ Mq (78) 



(in accordance with the result of AC) , where is give n in 
[85| and Npk is the average number density of peak ([54I Eq. 
4.4). Galaxies and Clusters of galaxies are correlated sys- 
tems whose autocorrelation function, Hr), can be expressed, 
in a power law form ([3i|;[8i];[83;@;[ii). The descrip- 
tion of dynamical friction in these systems need to use a 
distribution of the stochastic forces, W{F), taking account 
of correlations. In this last case the coefiicient of dynamical 
friction, 77, may be calculated using the equation: 



j d^FW{F)F^T{F)/{2 < 



>) 



(79) 



and using [83| distribution: 
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where A f , which is a linear integral function of the correla- 
tion function f(r), is given in the quoted paper (Eq. 36). 

1.9.3 Tidal torques and structure formation 

The explanation of galaxies spins gain through tidal 
torques was pioneered by Peebles (p^) performed the 
first detailed calculation of the acquisition of angular mo- 
mentum in the early stages of protogalactic evolution. More 
recent analytic computations ('o?!, [7ll, [HI; [93|; [ol]; [sl 
and numerical simulations (j96]) have re-investigated the 
role of tidal torques in originating galaxies angular momen- 
tum. One way to study the variation of angular momentum 
with radius in a galaxy is that followed by [58 ] . In this ap- 
proach the protogalaxy is divided into a series of mass shells 
and the torque on each mass shell is computed separately. 
The density profile of each proto-structure is approximated 
by the superposition of a spherical profile, i5(r), and a ran- 
dom CDM distribution, e(r), which provides the quadrupole 
moment of the protogalaxy. As shown by (ssj the net rms 
torque on a mass shell centered on the origin of internal 
radius r and thickness Sr is given by: 



(81) 



where qim, the multipole moments of the shell and a/m, the 
tidal moments, are given by: 

{q2mir)') = -^M!h I k^dkP (k) j2 {krf (82) 



(2^)^ 



(a2m(r)(72m('')> 



2plr- 



dkP [k] ji {krY 



(83) 



27r2' 



PbMsh / kdkP (k) ji (kr) j2{kr) 



(84) 

where Msh is the mass of the shell, ji (r) and j2 (r) arc the 
spherical Bessel function of first and second order while the 
power spectrum P{k) is given by, [s^l (equation (G3)): 



T{k) = 



[ln(l + 2.34g)] 
2.Mq 

■[l + 3.89q+ (m.lqY 



(5.46g)3 + (6.71)*]^5) 



Here 9 = pci7(l-686p^) represents 



(where q - Oxh^Mpc-i ■ """"^ " ~ a^"-/ v^-"""P7^ 
the ratio of the energy density in relativistic particles to 
that in photons {9 = 1 corresponds to photons and three 
flavors of relativistic neutrinos). The power spectrum was 
normalized to reproduce the observed abundance of rich 
cluster of galaxies. Filtering the spectrum on cluster scales, 
Rf = 3h~^Mpc, I have obtained the rms torque, r(r), on 
a mass shell using Eq. ([5T|) then I obtained the total spe- 
cific angular momentum, /i(r, u), acquired during expansion 
integrating the torque over time ([5g] Eq. 35): 



h{r, v) 



Tpto j-5/2 (1 - cos 



_ 1 

~ 3 V4 



(i9-sin</'/i(i?)-/2(i9)|s 



2/3 



d^ (86) 



the functions /i(i9), /2('?) are given by [58] (Eq. 31) while 
the mean over-density inside the shell, S{r), is given by [sst : 



where d{r) 



^^^^i^-^. As showed by [5^, the rms specific 
angular momentum, /i(r, v), increases with distance r while 
peaks of greater v acquire less angular momentum via tidal 
torques. This is the angular momentum-density anticorre- 
lation showed by [tH . This effect arises because the angular 
momentum is proportional to the gain at turn around time. 



tm, which in turn is proportional to (5(r, v) 2 
1.9.4 Modification of collapse 



-3/2 



Tidal torques and dynamical friction acts in a similar 
fashion. As previously reported, AC calculated the effect 
of dynamical friction taking into account the effect of sub- 
structure, showing that dynamical friction delays the col- 
lapse of low-j/ peaks inducing a bias of dynamical nature. 
Similarly non-radial motions would slow the rate of growth 
of the density contrast by lowering the peculiar velocity and 
suppress collapse once the system detaches from general ex- 
pansion. In fact, in the central regions of a density peak 
(r < 0.5Rf) the velocity dispersion attain nearly the same 
value while at larger radii (r > Rf) the radial component 
is lower than the tangential component. This means that 
motions in the outer regions are predominantly non-radial 
and in these regions the fate of the infalling material could 
be influenced by the amount of tangential velocity relative 
to the radial one. This can be shown writing the equation 
of motion of a spherically symmetric mass distribution with 
density n{r) ([971]): 

d d Til d 

-n{vr) + ^n{vf) + (2(«2) - (vD) -+n{r)-{vr) = 



where {vr) and (w^) are, respectively, the mean radial and 
tangential streaming velocity. Eq. ([551) shows that high tan- 
gential velocity dispersion ((w^) > 2(w^)) may alter the infall 
pattern. The expected delay in the collapse of a perturba- 
tion, due to non-radial motions, dynamical friction and also 
taking account of a non-zero cosmological constant, may be 
calculated solving the equation for the radial acceleration 

(HI; [11; AC; M)-- 



dvr L'^ir^v) 
~dt ~ MVS 



, , dr K , ^ 

9{r) - + Y (89) 



where L(r, v) is the angular momentum, g{r) the accel- 
eration, and A the cosmological constant. Writing the 
proper radius of a shell in terms of the expansion parameter, 
a(rj,t): 



remembering that 



r{r.i,t) = r,a{ri,t) 



47r 



M ^ ^Pbir,,t)a^r„t)rf 



(90) 



(91) 



and that pb = where Hq is the Hubble constant and 

assuming that no shell crossing occurs so that the total mass 
inside each shell remains constant, that is: 



a^(r^,t) 

the Eq. ([55)1 may be written as: 
d^a H^{l + 5) 4G2l2 da A 



dt^ 



2a2 



■n- 



Hi{l + SYrj"a^ dt 3 



(92) 



a (93) 
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FIG. 1: The time evolution of the expansion parameter. The 
solid line is a(t) for the SM; the dashed line is a(t) taking account 
only dynamical friction; the dotted line is a(t) taking account 
of the cumulative effect of non-radial motions and dynamical 
friction in the case oi & v = 2 peak. 



h{r,v) — L{r,v)/Msh found in Eq. ([86|l . to obtain a(t) 
and the time of collapse, Tc{r,v). As shown by [s^, this 
last quantity in the case of a pure SM (namely when tidal 
torques, dynamical friction and cosmological constant are 
not taken into account) is given by: 



Tcoir, v) 



(94) 



In Fig. 1, I show the effects of non-radial motions and dy- 
namical friction separately, in the case oi a. u = 2 peak. As 
displayed non-radial motions have a larger effect on the col- 
lapse delay with respect to dynamical friction. In Fig. 2, I 
compare the results for the time of collapse, Tc, for = 2, 
4 with the time of collapse of the classical SM (Eq. [M)) . 
As shown the presence of non-radial motions produces an 
increase in the time of collapse of a spherical shell. The 
collapse delay is larger for low value of v and becomes neg- 
ligible for V >'i. This result is in agreement with the angu- 
lar momentum-density anticorrelation effect: density peaks 
having low value of v acquire a larger angular momentum 
than high v peaks and consequently the collapse is more 
delayed with respect to high u peaks. 

Given Tc{r,v), I also calculated the total mass gravita- 
tionally bound to the final non-linear configuration. There 
are at least two criteria to establish the bound region to a 
perturbation 5{r): a statistical one ((59|). and a dynamical 
one ([131, summarized in the following. 

In biased galaxy formation theory structures form around 
the local maxima of the density field. Every density peak 
binds a mass M that can be calculated when we know the 
binding radius of the density peak. The radius of the bound 



^ As we told in the introduction, the inclusion of the peculiar velocity 
field changes the collapse as: HiTcO — ^ ^ \3/2 where c and ei 



2b 



FIG. 2: Fig. 2a. The time of collapse of a shell of matter in units 
of the age of the universe to for u = 2 (dotted line) compared 
with SM (solid line). Fig. 2b. The time of collapse of a shell of 
matter in units of the age of the universe to for v = A (dotted 
line) compared with SM (solid line). 



region for a chosen density profile 5{r) may be obtained in 
several ways. A first criterion is statistic. The binding ra- 
dius of the region, r^, is given by the solution of the equation: 



< 5{r) >^< {S- <d>f > 



1/2 



(95) 



([581]; [59|). At radius r << rt the motion of particles is 
predominant toward the peak while when r » rb the par- 
ticle is not bound to the peak. Another criterion that can 
be used is dynamical. It supposes that the binding radius 
is given by the condition that a shell collapse in a time, Tj,, 
smaller than the age of the universe tg: 



Tc{r) < to 



(96) 
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physics of the collapse process of a shell. For this reason I 
used it to calculate the binding radius. The time of collapse, 
Tc(r), at radius r can be obtained solving numerically Eq. 
([93]) for different values of Si, the initial overdensity, from a 
given density profile 6{r). I use the average density profile 
given by [53|: 



Sir) 



U(o)i/2 



7^(0)1/2(1-72) 



3 

(97) 



where A is a constant given by the normalization of the 
perturbation spectrum, P{k), v = where S^o = 1.686 

is the critical threshold for a SM, a{M) is the r.m.s. den- 
sity fluctuation on the mass scale Af, ^(r) is the correlation 
function of two points, 7 and R.^, two constants obtainable 
from the spectrum (see .52]) and finally 9{'yi>,j) is a func- 
tion given in the quoted paper (eq. 6.14). Given the average 
density profile the average density inside the radius r in a 
spherical perturbation is given by Eq. ((87| . 

Finally, I calculated the binding radius, rb{i'), for a SM, 
calculating Tco{r) by means of Eq. ([94]) and the density 
profile given in Eq. (|97p and then applying the condition 
Tcoir) < to- I repeated the calculation for 1.7 < v < 4. 
Then I repeated the calculation using T'c(r), the collapse 
time that takes into account non-radial motions and dy- 
namical friction. I found a relation between v and the mass 
of the cluster using the equation: M = ^r^pb. The result 
is the plot in Fig. 3 for the binding radius versus v. 

In f ig. 4, I compare the peak mass obtained from SM, 
using [571 criterion, with that obtained from the model tak- 
ing into account non-radial motions, dynamical friction and 
A ^ 0. As shown for high values of (:^ > 3) the two mod- 
els give the same result for the mass while for ly < 3 the 
effect of non-radial motions produces less bound mass with 
respect to SM. Decreasing the effect of non-radial motions 
produces a decrease in the bound mass. 

The situation represented in the previous three figures 
may be summarized as follows: dynamical friction and non- 
radial motions delays the collapse of perturbations. Both 
effects act in the direction of delaying structure collapse, so 
that their effects add. The effects have a similar magnitude, 
but non-radial motions induce a slight larger delay in col- 
lapse. As a consequence of this delay of collapse the matter 
bound to structures is less than what expected in the case 
of SM. 

1.9.5 The threshold of collapse Sc 

In this section, I am going to show how dynamical friction 
and tidal fields influence the critical overdensity threshold 
for the collapse, 5^ which is not constant as in a SM but it 
depends on mass. An analytic determination of Sdi') can 
be obtained following a technique similar to that used by 

Using Eq. (|93p it is possible to obtain the value of the 
expansion parameter of the turn around epoch, amax , which 



is characterized by the condition ^ = 0. 
between v and Si, in linear theory ([33|) 



Using the relation 
I find: 



B{M) 



rtaL2 ■ dr Aq 
GAPr^ l~fi{S) ' ^^6GM' 



where Sco — 1.686 is the critical threshold for SM, rj is the 
initial radius, rta is the turn-around radius, Xo = CqTco and 



FIG. 3: Variation of the binding radius rt with p. The solid 
line is the binding radius in the SM, while the dashed line is the 
same as in presence of non-radial motions dynamical friction and 
A / 0, for = 2. 



proto-structure during evolution that is calculated as shown 
in section 1.9.3. The result of the calculation is shown in Fig. 
5, where I plot Sdi') obtained by means of the model of the 
present paper together with that obtained by jlOOf (ST) us- 
ing an ellipsoidal collapse model. The dashed line represents 
Sc{iy) obtained with the present model, while the solid line 
that of ST. Both models show that the threshold for collapse 
decreases with mass and when v > S the threshold assume 
the typical value of the SM. In other words, this means that, 
in order to form structure, more massive peaks must cross 
a lower threshold, <5c(i^), with respect to under-dense ones. 
At the same time, since the probability to find high peaks is 
larger in more dense regions, this means that, statistically, 
in order to form structure, peaks in more dense regions may 
have a lower value of the threshold, Sdi'), with respect to 
those of under-dense regions. This is due to the fact that 
less massive objects are more influenced by external tides, 
and consequently they must be more overdense to collapse 
by a given time. In fact, the angular momentum acquired by 
a shell centered on a peak in the CDM density distribution 
is anti-correlated with density: high-density peaks ac quir e 
less angular momentum than low-density peaks ([7l|; |58|: 
[sgj). A larger amount of angular momentum acquired by 
low-density peaks (with respect to the high-density ones) 
implies that these peaks can more easily resist gravitational 
collapse and consequently it is more difficult for them to 
form structure. This is in agreement with [t^ . [63 |. which 
pointed out that the gravitational collapse is slowed down 
by the effect of the shear rather than fastened by it (as sus- 
tained by other authors). Therefore, on small scales, where 
the shear is statistically greater, structures need, on aver- 
age, a higher density contrast to collapse. This results in 
a tendency for less dense regions to accrete less mass, with 
respect to a classical SM, inducing a biasing of over-dense 
regions towards higher mass. 

1.10 Mass function 
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Finally trying to determine the luminos ity f unction obser- 
vatively is problematic (see for example |lOl| '). 
For the above reasons, the theoretical determination of the 
mass function is very important. One of the most success- 
ful study in the subject is that of [5l[ . This theory is based 
upon these hypotheses: 

• The linear density field is described by a stochastic 
Gaussian field. The statistics of the matter distribu- 
tion is Gaussian. 

• The evolution of density perturbations is that de- 
scribed by the linear theory. Structures form in those 
regions where the overdensity linearly evolved and 
filtered with a top-hat filter exceeds a threshold 5c 
(5^^ 1.68, obtained from the spherical collapse model 

• for (5 > 6c regions collapse to points. The probability 
that an object forms at a certain point is proportional 
to the probability that the point is in a region with 
S > Sc given by: 



FIG. 4: The mass accreted by a collapsed perturbation, in units 
of 10^^ Mq, taking into account non-radial motions, dynamical 
friction and a non-zero cosmological constant (dotted line) com- 
pared to SM mass (solid line). 




FIG. 5: The critical threshold, Scif) versus v. The dashed line 
is obtained with the model of the present paper while the solid 
hne is that of ST. 



by the relation: 



dN = n{M)dM 



(99) 



that is the number of objects per unit volume, having a mass 
in the range M ed M + dM. The multiplicity function can 
also be used to define the luminosity function after having 
fixed the ratio Obtaining the mass function starting 
from that of luminosity is complicated since the ratio 
is known with noteworthy uncertainty and it is different 
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The multiplicity function is given by: 
dP 

p{M, z) = -Pa-Q^dM = n{M)MdM 



^ (100) 



(101) 



If we add the conditions = 1, cx fc", the Press- 

Schechter solution is autosimilar and has the form: 



p{M,z) 
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where M^{z) cx {1 + z) 
the theory: 



Several are the problems of 



• Statistical problems: in the limit of vanishing 
smoothing radii, or of infinite variance, the fraction 
of collapsed mass, asymptotes to 1/2. This is a sig- 
nature of linear theory: only initially overdense re- 
gions, which constitute half of the mass, are able to 
collapse. Nonetheless, underdense regions can be in- 
cluded in larger overdense ones, or, more generally, 
non-collapsed regions have a finite probability of be- 
ing included in larger collapsed ones; this is commonly 
called cloud- in- cloud problem. PS argued that the 
missing mass would accrete on the formed structures, 
doubling their mass without changing the shape of the 
MF; however, they did not give a true demonstration 
of that. Then, they multiplied their MF by a "fudge 
factor" 2. Other authors used to multiply the MF by 
a factor (1 -I- /), with / denoting the fraction of mass 
accreted by the already formed structures. 

• Dynamical problems: the heuristic derivation of 
the PS MF bypasses all the complications related to 
the highly non-linear dynamics of gravitational col- 
lapse. Spherical collapse helps in determining the 6c 
parameter and in identifying collapsed structures with 
virialized halos. However, the PS procedure com- 
pletely ignores important dynamical elements, such 
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that every structure virializes just after collapse is a 
crude simplification: when a region collapses, all its 
substructure is supposed by PS to be erased at once, 
while in realistic cases the erasure of substructures 
is connected to the two-body interaction of already 
collapsed clumps, an important piece of gravitational 
dynamics which is completely missed by the PS pro- 
cedure. 

• Geometrical problems: to estimate the mass func- 
tion from the fraction of collapsed mass at a given 
scale it is necessary to relate the mass of the formed 
structure to the resolution 

In practice, the true geometry of the collapsed regions 
in the Lagrangian space (i.e. as mapped in the ini- 
tial configuration) can be quite complex, especially at 
intermediate and small masses; in this case a differ- 
ent and more sophisticate mass assignment ought to 
be developed, so that geometry is taken into account. 
For instance, if structures are supposed to form in the 
peaks of the initial field, a different and more geomet- 
rical way to count collapsed structures could be based 
on peak abundances. 

Despite all of its problems, the PS procedure proved suc- 
cessful, as compared to N-body simulations, and a good 
s tarti n g poi nt for all the subseque nt w o rks o n the su bjec t 
(lip; [103, Bond et al. (1991), [lo|, [lo|, [lO^, [lo3, 
[69l |. Most authors reported the PS formula to fit well their 
N-body results; nonetheless, all the authors agree in stat- 
ing the validity of the PS formula to be only statistical, i.e. 
the existence of the single halos is not well predicted by the 
linear overdensity criterion of PS (see in particular Bond et 
al. 1991)). Ther e are however some exceptions to this gen- 
eral agreement: jl08l | reported their MF, based on a CDM 
spectrum, to be very similar to a power-law with slope —2, 
different from the PS formula both at s mall and at large 
masses. Several authors, jl05| . jl09l | and |llC| noted that, 
to make the PS formula agree with their simulations (based 
on CDM or CHDM spectra), it is necessary to lower the 
value of the Sc parameter as redshift increases. The same 
thing was found by lll| , but was interpreted as an artifact 
of their clump-finding algorithm. Recent simulations seem 
to confirm this tr end. 

Lacey & Cole ([Til]) extended the comparison to N-body 
simulations to the predictions for merging histories of dark- 
matter halos; they found again a good agreement between 
theory and simulations. This fact is noteworthy, as merging 
histories contain much more detailed information about hi- 
erarchical collapse. Several improvements of the theory ex- 
ists: [Toil (the so called Extended Press & Schechter (EPS) 
formalism; jll3l |). 

1.11 CDM, HDM and others cosmogonies 

The study of origin and formation of structure in the uni- 
verse has been historically fundamentally framed into two 
theories: the CDM theory, in which WIMPS constitutes 
the main part of Dark Matter, and HDM in which neutri- 
nos dominate. As we are going to see, structure formation 
in these scenarios is completely different since WIMPS and 
neutrinos are subject to different physical phenomena and 
then the transfer function is noteworthy different in the two 
cases. Both theories have the same starting points: 
1) The universe is fundamentally constituted by Dark Mat- 
ter (WIMPS in the CDM model and neutrinos in the HDM) 
and p — Pc { ft — 1). 



3) Fluctuations originating in the primordial universe are 
adiabatic, scale invariant, n = 1, and Gaussian ((28|). 
If universe is dominated by neutrinos with mass m,y = 30ei; 
the transfer function is determined by the free-streaming (or 
Landau damping) of neutrinos. This phenomenon consists 
in the smoothing of inhomogeneities in the primordial uni- 
verse (due to perturbations in the acoUisional components) 
because of the motion of neutrinos from overdense to un- 
derdense regions. Neutrinos diffusion and the smoothing 
of inhomogeneities is possible only before t = teq- After 
this epoch, there is no longer free-streaming but the density 
perturbation has definitely changed by its previous action. 
Free-streaming scale or mass can be estimated calculating 
the distance covered from a particle decoupled from plasma. 
Results that one obtains for the free-streaming scale and 
mass is ([Tsj): 



A 



FS 



mMpc{m^/mevy 



M 



FS 



lO^%m^/30evy^Mo 
(103) 

Because of free-streaming of neutrinos the HDM spectrum 
is characterized by a cut-off at short wave-length. The re- 
sult is that the first objects that form are superclusters and 
structure formation proceeds because of fragmentation. [47j 
showed that the first structure to form are flat and were 
called "pancakes". After formation, these objects enter in 
the non-linear phase along one of the axes and baryons in- 
side start to collide and dissipate their gravitational energy. 
Galaxies form for fragmentation processes. Structure for- 
mation follows a 'top-down' scheme, that means that larger 
objects (e.g., clusters) form before, while smaller objects 
( e.g., g alaxi es) later. N-body simulations of HDM universes 
f |ll4| : jll5l | showed that on scales larger than 10 Mpc struc- 
tures is qualitatively similar to voids and to the filamentary 
structure that is visible in the CFA, but the clustering mea- 
sured in N-body simulations is larger than that observed in 
the CFA. When one tries to reproduce the observed correla- 
tion function, one arrives to the conclusion that pancaking 
should have happened at redshift z < 1, in disagreement 
with observed galaxies having z > 1 and QSO with z > 3. 
A further problem of the model is that of the peculiar veloc- 
ity that are smaller than values obtained from observations 
(111)- 

The HDM model after a series of studies in the '80s has been 
abandoned for the problems it has and replaced by another 
model, the CDM, which is in better agreement with obser- 
vations. The CDM model has a spectrum without a cut-off 
at short wave -length (at least till scales much smaller than 
galactic scales) because the damping scale is unimportant 
for WIMPS with mass > 1 Gev. Structure formation is 
typically hyerachical: from smaller scale structure to larger 
ones. This scheme is a 'bottom-up' scheme. When the CDM 
model was introduced it obtained noteworthy successes in 
the description of the characteristic of the universe (cluster- 
ing statistics of galaxies, peculiar v elocit ie s, C MBR fl uctu - 
ation s) from the galactic scale on ([116*1; [ll7l |: (s^ : jll8l |: 
|ll9l |: |28j). The model has shown some weak points, when 
compared with more and more precise data. 
The reason of the success of the CDM model is fundamen- 
tally due to the fact that WIMPS interact with matter by 
means of gravity only, and does not feel the effect of pressure 
forces due to interaction with radiation (to which matter 
components are subject). Structure formation starts before 
in the CDM component, at t < teq, which give rise to the 
potential wells in which baryonic matter can then fall. It 
is important to notice that, in order to reproduce observa- 
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mass. ^. Typical problems of the CDM model in absence 
of biasing are the too high values for the r.m.s. of pecu- 
liar velocity of couple of galaxies (values of 1000 km/s vs. 
300 ± 50 km/s observed). Another problem is the correla- 
tion length, ro, in N-body simulations, for the correlation 
function ^ which is equal to 1.3h~^Mpc, smaller than the 
observed value 5.5h~^Mpc. On the other side, if, in order 
to eliminate the quoted problems one introduce the biasing 
hypothesis, there is the supplementary problem of finding a 
physical mechanism that explai ns th e origin of bi as. Several 



121|; [123) but there 



conjectures has been proposed ( |12C 
is no full agreement on them. 

Summarizing, one can tell that although at the beginning 
the standard form of CDM was very successful in describ- 
ing the structures observed in the Universe (galaxy cluster- 
ing statistics, structure formati on ep ochs, peculiar velocity 
flows) ([ii^; [113; [El; [ni; [ui; HI) recent measure- 
ments have shown several deficiencies in the model, at least, 
when any bias of the distribution of ga la xies r elativ e to th e 
mass is constant with scale (see jl23j |: |l24| : jl25j . [l26j |). 
Some of the most difficult problems that must be reconciled 
with the theory are: 

• the magnitude of the dipole of the a ngular distribution 
of optically selected galaxies f jl27l |): 

• the possible observations of cluster s of galaxies with 
high velocity dispersion at z > 0.5 ( [l28| |): 

• the strong clustering of rich clusters of galaxies, 
^cc(r-) ~ (r/25/i"^Mpc)"^, far in excess of CDM pre- 
dictions (|8y]); 

• the X-ray temperature distribution function of clus- 
t ers^ over-producing the observed clusters abundances 

• the conflict between the normalization of the spectrum 
of the perturbation which is required by different types 
of observations; 

• the incorrect scale dependence of the galaxy correla- 
tion function, ^(r), on scales 10 to 100 h~^Mpc, hav- 
ing ^(r) too little power on the lar ge sc al es co m pared 
to th e power on smaller scales ( |l29l |: [l30j |: [l3lj |: 



• Normalization obtained from COBE data ( |133j ) on 
scales of the order of 10'^ Mpc requires as = 0.95 ±0.2, 
where erg is the rms value of in a sphere of 
8h~^Mpc. Normalization on scales 10 to b OMpc ob- 
tained from QDOT and POTENT ([HI) requires 
that CTg is in the range 0.7 1.1, which is compati- 
ble with COBE normalization while the observations 
of the pairwise velocity dispersion of galaxies on scales 
r < SMpc seem to require erg < 0.5. 

• Another problem of CDM model is the incorrect 
scale dependence of the galaxy correlation function, 
^(r), on scales 10 to 100 Mpc, having ^(r) too little 
power on the large scales compared to th e power 
on smaller scales. The APM survey ( [l29j |'). giving 
the galaxy angular correlation function, the 1.2 J y 
IRAS power spectrui n, t he QDOT survey f [l30j |). 
X-ray observations f [l35l |) and radio observations 



([13l|; [132]) agree with the quoted conclusion. As 
shown i n rec e nt stu dies of galaxy clustering on large 
scales ( [l29l |: [l30l |) the measured rms fluctuations 
within spheres of radius 20h~^Mpc have value 2- 
3 times larger than that predicted by the CDM model. 



• Density profiles of CDM halos: the cusp obtained from 
numerical simulations seems too steep. - 

• simulations might yield too many satellites for galax- 
ies like our own. Though this second problem may 
have been the result of bad comparison of simulations 
with observations. This yielded a surge of interest in 
the last 2-3 years for Warm Dark matter and slightly 
collisional matter. 

These discrepancies between the theoretical predictions of 
the CDM model and the observations led many authors to 
conclude that the shape of the CDM spe ctrum is in c orrec t 
and toseaxch alternative models ([136[; [l37j : [l38l |: [l39l |: 

i; EH). 

Alternative models with more large-scale power than CDM 
have been introduced in order to solve the latter problem. 
Several authors have lowered the matter density under the 
critical value {flm < 1) {^m — 0.3) (this model is called 
open cold dark matter model (OCDM)) and others f j l36l |: 
Efstathiou et al. 1990a; ^3,]) have also added a cosmological 
constant in order to retain a flat Universe (f2„i + ~ 
This model is known as ACDM model. The spectrum of 
the matter density is specified by the transfer function, but 
its shape is affected because of the fact that the epoch of 
matter-radiation equality is earlier, 1 -I- Zeq being increased 
by a factor 1/flm- The epoch of matter-radiation equality 
is earlier, because 1 -I- z^q is increased by a factor l/rim. 
Around the epoch z\ the growth of the density contrast 
slows down and ceases after z\. As a consequence the 
normalization of the transfer function begins to fall, even 
if its shape is retained (and pushes its impri nt to l arge r 
scales). Mixed Dark Matter models (MDM) ([l33|; [T3a |: 
[l39l |: [l40l |) increase the large-scale power because neutrinos 
free-streaming damps the power on small scales. Alterna- 
tively, changing t he pr imeval spectrum several problems of 
CDM are solved ([mH). For example, in the rCDM model 
the needed changes in the power spectrum may be obtained 
in = 1 CDM models if matter-radiation equivalence is 
delayed, such as by the addition of an additional relativis- 
tic particle species. Finally it is possible to assume that 
the threshold for galaxy formation is not spatially invariant 
but weakly modulated (2% — 3% on scales r > lOh^^Mpc) 
by large scale density fluctuations, with the result that the 
clustering on large-scale is significantly increased ( jl24{ '). 

1.12 Constraints from recent astrophysical observations 

In the past decade we have witnessed spectacular progress 
in precision measurements in astrophysics as a result of sig- 
nificant improvements in terrestrial and extraterrestrial in- 
strumentation. The (second phase of the) Hubble telescope 
opened up novel paths in our quest for understanding the 
Universe, by allowing observations on distant corners of the 
observable Universe that were not accessible before. 

From the point of view of interest to particle physics, the 
most spectacular claims from astrophysics came 8 years ago 
from the study of distan t sup ernovae (redshifts z ~ 1) by 
two independent groups [l42l |. These observations pointed 
towards a current era acceleration of our Universe, some- 
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Universe, or in general by a non-zero dark energy compo- 
nent, which could even be relaxing to zero (the data are 
consistent with this possibility). This claim, if true, could 
revolutionize our understanding of the basic physics gov- 
erning fundamental interactions in Nature. Indeed, only 
a few years ago, particle theorists were trying to identify 
(alas in vain!) an exact symmetry of nature that could set 
the cosmological constant (or more generally the vacuum 
energy) to zero. Now, astrophysical observations point to 
the contrary. The skeptics may question the accuracy of 
the supernovae observations, however, there is additional 
evidence from quite different in origin astrophysical obser- 
vations, those associated with the measurement of the cos- 
mic microwave background radiation (CMB), which point 
towards the fact that 73 % of the Universe vacuum energy 
consists of a dark (unknown) energy substance, in agreement 
with the (preliminary) supernovae observations. Moreover, 
recently [14S | two more distant supernovae have been dis- 
covered (z > 1), exhibiting similar features as the previous 
measurements, thereby supporting the geometric interpre- 
tation on the acceleration of the Universe today, and arguing 
against the nuclear physics or intergalactic dust effects. 

Above all, however, there are the very recent data from 
a new probe of Cosmic Microwave Background Radia- 
tion Anisot ropy (Wilkinson Microwave Anisotropy Probe 
(WMAP)) [144| |. In its first year of running WMAP mea- 
sured CMB anisotropics to an unprecedented accuracy of 
billionth of a Kelvin degree, thereby correcting previous 
measure ment s by the Cosmic Background Explorer (COBE) 



satellite |145| by several orders of magnitude. This new 
satellite experiment, therefore, opened up a new era for 
astroparticle physics, given that such accuracies allow for 
a determination (using b est f it models of the Universe) of 
cosmological parameters [l46j . and in particular cosmolog- 
ical densities, which, as we shall discuss in this review, is 
quite relevant for constraining models of particle physics to 
a significant degree. 

The WMAP satellite experiment determined the most im- 
portant cosmological paramete rs th at could be of relevance 
to particle physicists, namely [l46j |: the Hubble constant, 
and thus the age of the Universe, the thickness of the last 
scattering surface, the dark energy and dark matter content 
of the Universe (to an unprecedented accuracy) , con firm- 
ing the earlier claims from supernovae la data jl42l |. and 
provided evidence for early reionization (z ~ 20), which, at 
least from the point of view of large scale structure forma- 
tion, excludes Warm Dark Matter particle theory models. 

An important comment concerns the dark energy com- 
ponent (73 % ) of the Universe. The WMAP measured 
equation of state for the Universe p = wp, with p the pres- 
sure and p the energy density, implies — 1 < w < —0.78 
(assuming the lower bound for the oreti cal reasons, other- 
wise the upper limit may be larger [l46j ). For comparison 
we note that w = —1 characterizes a perfect fluid Universe 
with non-zero, positive, cosmological constant. As we shall 
remark, supergravity quintessence models do have this fea- 
ture of w —1, and it may well be that by exploiting 
further the data on this dark energy component of the Uni- 
verse one may arrive at the physically correct supergravity 
model which could constrain the supersymmetric particle 
physics models. 

The results of the WMAP analysis (alone), including di- 
rectly measurable and derived quantities, are summarized 
in the tables appearing in figures I7I8I 

One therefore obtains the chart for the energy and matter 
content of our Universe depicted in figure 6. This chart is 




FIG. 6: The energy content of our Universe as obtained by 
fitting data of WMAP satellite. The chart is in perfect agreement 
with earlier claims made by direct measurements of a current 
era acceleration of the Universe from distant supernovae type la 
(courtesy of http : //map .gsf c . nasa. gov /) . 



Table 1 

Power Law ACDM Model Parameters- WMAP Data Only 
Parameter Mean (68% confidence range) Maximum Lil<eliliood 



Baryon Density 




0.024 ±0.001 


0.023 


Matter Density 




0.14 ±0.02 


0.15 


Hubble Constant 


h 


0.72 ± 0.05 


0.68 


Amplitude 


A 


0.9 ±0.1 


0.80 


Optical Depth 


T 


0.166iro?; 


0.11 


Spectral Index 




0.99 ±0.04 


0.97 

1431/1342 



"Fit to WMAP data only 



FIG. 7; Cosmological parameters me asured by WMAP 
(only):directly measurable quantities |l46l ]. 



novae la Data [142 |. 

It should be stressed that the interpretation of the su- 
pernovae data is based o n a best fit Friedmann-Robertson- 
Walker (FRW) Universe [m^: 

- o.er^A ^ -0.2 ± 0.1 , for vlm < 1.5 (104) 

with flM,A corresponding to the matter and cosmological 
matter densities. Assuming a flat model (k=0), fitotai = Ij 
as supported by the CMB data, the SNIa data alone imply: 



* = 0.28^°;°^ {la stat.)l'^-^^l {identified syst.) (105) 



+0.05 



^ , where a 



The deceleration parameter defined as q e 
is the cosmic scale factor, receives the following form if we 
omit the contribution of photons which is very small, 



-0.57 < 



{^A 0.7) . (106) 



Hence (jl04p and (|105p provide evidence for a current era 
acceleration of the Universe. At this stage it should be 
stressed that the recent observation of two more supernovae 
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Table 1 

Derived Cosmological Parameters 



Parameter 


Mean (68% confidence range) 


Amplitude of Galaxy Fluctuations 


(Ts = 0.9 ±0.1 


Characteristic Amplitude of Velocity Fluctuations 


asfllf = 0. 44 ± 0.10 


Baryon Density/Critical Density 


116 = 0.047 ±0.006 


Matter Density/Critical Density 


n„, =0.29 ±0.07 


Age of the Universe 


to = 13.4 ±0.3 Gyr 


Redshift of Reionization (ionization fraction, Xe = 1) 


:^ = 17±5 


Redshift at Decoupling 


Zdec = IO88I2 


Age of the Universe at Decoupling 


tj^^ = 372 ± 14 kyr 


Thickness of Surface of Last Scatter 


Azdcc = 194 ± 2 


Thickness of Surface of Last Scatter 


Atd^^ = 115 ± 5 kyr 


Redshift at Matter/Radiation Equality 




Sound Horizon at Decoupling 


= 144±4 Mpc 


Angular Diameter Distance to the Decoupling Surface 


dA = 13.7 ±0.5 Gpc 


Acoustic Angular Scale(/_4 = Trd^/r^) 


a A = 299 ± 2 


Current Density of Baryons 


Tib = (2.7 ± 0.1) X lO""^ cm"^ 


Baryon/Photon Ratio 


■r, = (G.rjtoi) X 10-1" 



FIG. 8: Cosmological par ameters measured by WMAP 
(only) :derived quantities [l46l ]. 

and argues rather against the role of nuclear (evolution) or 
intergalactic dust effects. 

The recent data of WMAP satellite lead to a new determi- 
nation of Sltotai = 1.02 ± 0.02, where ^itotai = Ptotai/Pc due 
to high precision measurements of secondary (two more) 
acoustic peaks as compared with previous CMB measure- 
ments. Essentially the value of f2 is determined by the posi- 
tion of the first acoustic peak in a Gaussian model, whose re- 
liability increases significantly by the discovery of secoi idary 
peaks and their excellent fit with the Gaussian model [l46j . 

2. CONCLUSIONS 

This paper provides a review of the variants of dark mat- 
ter which are thought to be fundamental components of the 
universe and their role in origin and evolution of structures. 
It moreover gives some new original results concerning im- 
provements to the spherical collapse model. In particular, 
how the spherical collapse model is modified when we take 
into account dynamical friction and tidal torques. Studies of 
several decades have shown that, if we have a right knowl- 
edge of the law of gravity, dark matter is a fundamental 
component of our universe. While models based upon Hot 
Dark Matter (e.g., neutrinos) gives a reasonable descrip- 
tion of structures on large scales models based upon Cold 



Dark Matter (e.g., axions) are more successful in describing 
small and intermediate scales. A fundamental ingredient in 
the recipe of structure formation is inflation which provides 
a spectrum of adiabatic Gaussian perturbations which can 
be well described by a power-law spectrum, tilted from the 
Harrison-Zel'dovich spectrum, normally tilted so as to pro- 
vide extra large scale power. The magnitude of the tilt may 
be modest or pronounced. The details of structure forma- 
tion are very sensitive to the matter content of the universe. 
It appears that if cold dark matter is the main constituent 
of the universe, present observations require that the initial 
perturbations be adiabatic — isocurvature perturbations 
generate excessively large cmb anisotropics for the same fi- 
nal density perturbation. Adiabatic perturbations are ex- 
actly what infiation provides. In CDM models, the only 
remaining alternative would appear to be texture seeded 
models, which have been placed in jeopardy by a combina- 
tion of microwave anisotropy and velocity data, though the 
death blow apparently remains to be struck. The survey was 
completed by examining variants on the CDM model which 
may be better suited to explaining the observational data. 
The standard technique is to utilize additional matter (be 
it a component of hot dark matter or of a cosmological con- 
stant) to remove short-scale power from the CDM spectrum. 
Hot dark matter does this by free-streaming, a cosmologi- 
cal constant by delaying matter-radiation equality. Because 
this power can be removed over a much shorter range of 
scales than with tilt, it enables an explanation of the ob- 
served deficit of short-scale power relative to intermediate 
scale power in the spectrum. 

MDM (Mixed dark Model) adds yet another new parame- 
ter, roughly speaking an ability to remove short-scale power 
from the spectrum while leaving large scales untouched, and 
may be necessary should all present observations stand up. 
It appears likely that MDM will however need an initial 
spectrum close to n = 1 wit h no gravitational waves if it 
is to succeed. Some studies ( |l36l |: Efstathiou et al. 1990a; 
[^) has shown that lowering the matter density under the 
critical value (f2m < 1) and adding a cosmological constant 
in order to retain a flat Universe {ilm + = 1), gives good 
results in the c ase o f flm — 0.3. Moreover new observational 
evidences (see |l47( ) indicates that we are living in a A 7^ 
universe. 
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Table 3: Estimates of i J,„ 



Global Measures 



Intlation. Occam 
Lmnimositv distance SNla 



Lens Counts 
CMS Peat 



" See Dekel et al. 1996 for a discussion of the table, 
is the rms mass density fluctuation in a top-hat sphere of 
radius ^h-'^Mpc. 

^ ,3 = S 6/ for IRAS galaxies, 6d for optical galaxies. 

' K/b, ~ 1.3, i)„ ~ l/ffg. 

il,a + !iA = 1 = 1, !!a = 0) 

-0.3 < f!,„ - S"!a < 2.5 (90%) 

Flat S!,„ > 0.19 (95%) 

Flat S!,„ > 0.3 1 (95%,) 

!!„. + < 1.5 (95%) 

!!„ + !!a > U.3 (likely ~ 0.7) 

S!„ - O.TSiA < 1.3 (likely < 0) 



Vh-ialized Objects (MIL)C 

Baryon fraction 
Cosmic Virial Th. 
Local Group 

Large-Scale; Flows Peculiar velocities 

Redshift Distortions 
Velocity vs Density 

Cluster Abundance 



a,a ~ 0.25 (0.1 - 1.0) 
Simfief ~ 0-3 - 0.5 (low-high Slj) 
Point mass £lm ~ 0.2 (halos — > 1) 
Pomt mass fim~0.15 (halos^O.7) 

!!„ > 0.3 (95%) 

Si^'a? = 0.8 ± 0.2 (/3J ~ 1.05«) 
A -0.5 -1.2 

A— 0.5 — 1.2 (scale dependent) 
/3„~ 0.4 -0.95 

iilfas = 0.5 - 0.6 (ft ~ 0.7 - 0.8=) 



Fluctuation Growth Cluster Morphology 
Galaxy Formation 
ft(rt vs C, 
Pi,(v) vs C, 



!!„ > 0.2 (?) 
(?) 

CDM n = 1 5 = 1: !!„/i66 ~ 0.3 
CDM flat: !J„/i6sn^ ~ 0.7 ± 0.1 




2 4 6 8 10 

r (Mpc) 



